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Digital Computer Laboratory 
Massachusetts Institute of Technology- 
Cambridge 39 t Massachusetts 



SUBJECT: BIWEEKLY REPORT. MARCH 21. 195** 

To: Jay W. Forrester 

From: Scientific and Engineering Computation Group 

1. MATHEMATICS. CODING AND APPLICATIONS 

l«i Introduction 

During the period covered by this report 33^ coded programs were run on 
the time allocated to the Scientific and Engineering Computation (S&EC) 
Group. These programs represent part of tie work that has been carried on 
in 32 of the problems that have been accepted by the S&EC Group* Progress 
on 20 of these problems is given below in terms of programming hours, minutes 
of computer time, and progress reports as submitted by the programmers in 
question* 

Two new problems were initiated during this period. Problem #1?1 will 
investigate modifications to existing techniques for the numerical calcul- 
ation of Fourier cosine transforms to obtain improved power spectrum esti- 
mates. This investigation is being carried out by D.T.Ross in partial 
fulfillment of the requirements for an M.S» degree in the Electrical Eng- 
ineering Department. In problem #1?2, F.J.Corbato' of the MIT Physics 
Department is developing routines for the evaluation of two-center overlap 
integrals that arise in molecular and crystal physics. 

H.B.Frost has completed the calculations he needed for his study (#152) 
of diffusion in an oxide coated cathode. A. summary of this problem will 
appear in Whirlwind Summary Report No. 37 • ^he results will be included in 
Mr. Frost's doctoral thesis to be submitted to the MIT Electrical Engineering 
Department. 

Intensive work has continued on the revision of the Comprehensive System 
of Service Routines (CS II) ♦ A brief internal memo (Bulletin Board Memo §55) 
has been circulated defining CS II and indicating how it differs from CS I. 
A detailed description of these changes will soon be available in Digital 
Computer Laboratory Memorandum M-2741. Use of CS II has already begun. How- 
ever, it is not expected to be used on a routine basis until after April 1, 195^« 
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1.2 Programs and Computer Operation 

The following summary is included as a guide for interpreting the 
abbreviations used below. A more detailed description of the terms involved 
can be found in M-2497. 

a. The upper case letter following the problem number has the following 
significance: 

A implies the problem is NOT for academic credit, is UN sponsored. 
B implies the problem is for academic credit, is UNsponsored. 
C implies the problem is NOT for academic credit, IS. sponsored. 
D implies the problem is for academic credit, IS sponsored. 

The absence of a letter indicates that it is an internal S&EC problem, 

b. DIC denotes the Division of Industrial Cooperation. 
DCL denotes the Digital Computer Laboratory. 

CMMC denotes the Committee on Machine Methods of Computation. 
DDL denotes the Division of Defense Laboratories. 

100. Comprehensive System of Service Routines, developed by the S&EG 

Group at the Digital Computer Laboratory for the input conversion 
of suitably prepared punched paper tapes. When so requested , 
these routines automatically provide a program with suitable 
programmed arithmetic, cycle-counting, and output facilities. 
sDGL Staff: Arden, 26 hours; Gombelic, 65 hours; Demur jian, 4-2.25 
hours; Denman, 82 hours; Frankovich, 45 hours; Helwig, 40 hours; 
Kopley, 8 hours; Porter, 18 hours; Siegel, 56 hours; Wfl, 1561 : 
minutes 

Testing and modification of the expanded system of utility programs 
has been continuing. 

The input program, the summer session program and the various test 
programs appear to be working correctly. 

The CS II programs are being modified so that the programs in the 
CS are called in from the drum instead of magnetic tape unit 0. The programs 
appear to work correctly up to the third pass. The output adaptation programs 
and the output blocks will be added after the program has been successfully 
tested. 



Testing of the post-mortem program continues. 



Staff 



The automatic post-mortem program for CS II has been undergoing 
test. Several minor errors in the program have been detected and will be 
corrected. 

Some difficulty was encountered with the buffer drum, which was not 
operating properly. This drum malfunction caused several correctly written 
programs to fail to work and resulted in the loss of considerable computer 
time. 

Siegel 
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Generalized decimal numbers have been included as a possible auto- 
matic post-mortem mode and a programmed arithmetic post-mortem is being 
tested in preparation for inclusion as an automatic feature. Testing of 
these and other minor changes is continuing, 

Arden 

A new PA Post-Mortem (PA PM) program has been written for CS II,. It 
forms a part of the new system of post-mortem programs to be stored on Magnetic 
Tape Unit and semi-permanently on the buffer drum. In response to a request 
for a PA M on the Manual Intervention Register (MIR), this PA PM program 
examines the stopped program (which has been stored on group of the auxiliary 
drum) for the presence of a PA. If it is a CSI PA, the CS I PA PM is given 
(the CS I PA PM program is being rewritten for this purpose by D. Combelic), 
and control returns to the main HI control program. If neither a CS I nor 
CS II PA appears to be in the stopped program, *'no PA" is printed out on the 
indicated output unit. If a CS II PA is found, the program continues with 
the CS II PA PM, 

If this is a tape PM request, the program again looks for a PA, and 
if it finds a CSI PA, gives this PA HI, If neither PA is found, the program 
i searches the requests appearing on the tape, and if either interpreted instruct- 
ions or generalized decimal numbers have been requested somewhere on the 
request tape, M no PA n is again printed out, -If neither of these types of 
requests is found, nothing is printed. In either case, the program then 
continues to give the PA f s requested on the tape. 

If a PA PM is requested via the MIR or is initiated by a tape request, 
and a CS II PA is found, but has not been used (the CS II PA is being modified 
so that it can give this information), the words M PA unused" are printed, and 
control goes back to the main PM control program. 

If a CS II PA is present and used, the PA PM program prints out the 
number system used in the PA, a title, where the program stopped, what instruct- 
ion was being performed when the program stopped, the contents of the register 
referred to by this instruction ( if this is significant), the contents of 
the MRA, the absolute address of the first register of the first buffer used 
and the numbers and contents of all the buffers (if no buffers are used, this 
is omitted) , similarly for the cycle-counters , and the jump table (the last 
•5 transfers of control by icp or isp instructions). 

A typical CSII PA PM print might appear as follows: 

(24 6) PA PM 

stopped at 105 105lits300 300 1 +.12345678 1-06 MRa| +.1234567891 -23 

l?65lfol + ol23456789|-12 lb I -.123456789 1 +37 

177l|o|o,5 l|2,7 2lo,0 3^,11 4l!|l*,3 5\ 2,6 

57lisp64 66licp4l 52licp54 57lisp64 ' 77licp?2 

* The. three vertical bars( HI) after a counter indicate that this was the last 
counter referred to by the program. 

Denman 
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The output blocks for CS II were rewritten to permit changing the 
location in storage of the block length table. The original location would 
have limited the number of output requests that could be processed by the 
output adaptation . 

In CS II the request tt MOAV tt will record a decimal point on the delayed 
printer. The stop code, for which this request was used in CS I, will be 
handled by two new special requests. They are n KQk end" for regular Wl mode 
and "iMOA end" when exit is desired in the interpreted mode. "MQA. end" and 
"iMOA end" record successively on magnetic tape a shift to lower case, stop 
code, two carriage returns and another stop code character-. These requests 
should be used to mark the end of recording on magnetic tape output. 

Demur jian 

Step 3 of the Output Adaptation for CS II has been translated into 
an optimum-length program by removing provisions which are no longer required. 

Kopley 

The final two tests for the first step of the Output Adaptation were 
run successfully. This program is now ready for. inclusion in CS II. 

The fourth step of the Output Adaptation has been rewritten and is 
now ready to be tested. This step has been revised so that it will work 
with either direct or indirect conversion. 

Porter 

101 C. Optical Properties of Thin Metal Films on transparent backings 

are determined and printed out, automatically by v this program; 
the input data consist of the observed reflection and trans- 
mission coefficients, the index of the backing, the wavelength, 
and the sample thickness. The program calculates by means of 
an iterative procedure and prints out the index of refraction 
and the absorption coefficient of the film, the rate of var- 
iation of these constants with reflection and transmission, 
and the film's conductivity and dielectric constant, 
s for Professor L. Harris, Chemistry Department, Dr. A.L.Loeb 
sby Dr. A..L..Loeb(DlC), 15 hours; Richmond; 30 hours 
:DCL Staff: Denman, 1 hourf HWI, 23 minutes 

In the last biweekly report equations were mentioned which enable 
one to obtain a first estimate of the optical and electrical constants of 
thin metal films for infrared radiation. These equations have been com- 
bined with the program computing reflection, and transmission in terms of 
the optical constants to check the validity of the approximate equations. 
Thus a set of R and T is ,taken to begin with, and the optical constants 
are calculated approximately. These constants are then stored in the 
"Main Program" for calculation of the corresponding R and T. If the differ- 
ence between the calculated and observed values is less than the expected 
experimental error, the approximation is considered satisfactory. The 
entire set of data reported by Waltersdorff has been run in this manner 
in about six minutes, and the results showed that the approximate equations 
are good to excellent. 
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106 C. MIT Seismic Project is concerned ( With the development of methods 
for locating deep reflections from underground, strata in 
seismic prospecting. The basic method is one of prediction 
by means of an optimum linear operator* 
:for Professor PJl.Hurley, Geology 'and Geophysics; Professor 
G. Wads-worth,, Mathematics Department 

: by E .A .Robins on( Res, Assoc .); Briscoe, J6 hours; Simpson, 20 hours; 
Walsh, 25 hours 
:DGLs WWI, 357 minutes 

Recent work on problem 106 has been concentrated on the computation 
of frequency information from seismograms, in the form of either power spectra 
or Fourier amplitudes and phases. This |nf ormation is necessary for the 
consideration of certain new approaches to linear operator selection. In 
one program a density plot scope display is used as an output ♦ 

108 C. An Interpretive Program is being developed that will accept algebraic 
equations, differential equations, etc. expressed on Flexowriter 
punched paper tape in ordinary mathematical notation (within 
certain limits imposed by the Flexowriter) as input and auto- 
matically provide the desired solution. 
:for Dr. J. H.Laning, Jr., Instrumentation Laboratory 
:by Dr. J.H.Laning, Jr. (DIC) 4 10 hours; Zierler, 40 hours; Block, 45 
hours 
:DCL : WSI, 156 minutes 

The runs during the^ past two weeks have been primarily testing 
problems. The function routines F^-F^" (see Laning and Zierler^ Instrument- 
ation Laboratory Report E-364) have been tested and work. An attempt to 
use the Photoelectric Tape Reader instead of the mechanifcall tape reader 
to read in the problem tape (line by line) has not been successful. 

A problem for T.X.Toong (see #-2699 Biweekly - February 21, 1954), • 
after a programmer's error hads&een corrected, ran successfully giving 
very satisfactory data. 

113 C. A Stress Analysis of an L-shaned Homogeneous Planar Structure is 
being made for the case of a concentrated static load. This 
structure is approximated by a framework of bars which will 
deform in the same manner as the prototype. This framework 
is then analyzed using the principles' of virtual work and 
Soutbwell relaxation techniques * Boundary conditions have 
been specified for the edge of the framework so that the 
deformations of the model will conform to the actual de- 
r formations of the structure. 

:for Professor J. S •Archer, Department of Civil and Sanitary Eng. 

:by S. Sydney(Res. Assist p GMMC), 35 hours 

:DCL: WWI, 186 minutes 

An analysis has been made for two load conditions on the "coarse" 
grid. Stress patterns have been obtained for these two cases, and these 
results will be inserted into the "fine" grid as initial conditions. This 
technique taay reduce ■ ^ the problem of convergence since the initial 
data will be closer to the final results* 
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119 C. Spherical Wave Propagation produced by the sudden release of a 

spherical distribution of compressed air in the atmosphere is 
being studied by numerical means, This involves replacing a 
set of non-linear hyperbolic partial differential equations 
in 2 independent arid 2 dependent variables by a set of differ- 
ence equations written along characteristics. An iterative 
procedure is used, to solve these equations* 

:for Professor C.C.Lin, Mathematics Department 

:by A. Ralston(CMMC), 15* hours 

:DCL: WWI, 252 minutes 

Further trouble in the past two weeks was caused by the noneonver- 
gence of the numerical process at one point. This caused the program to • 
get into a loop. This trouble has now been corrected and a production run 
is scheduled for early in the next biweekly period. 

120 D. The Aerothermopress or . This problem is concerned with the devel- 

opment of a device for increasing the stagnation pressure of 
a hot^high velocity gas stream by means of evaporative cool- 
ing, the analytical investigation being carried out on WWI 
involves the step-by-step solution ! of seven simultaneous non- 
linear differential equations which describe the thermodynamic 
and dynamic behavior of the compressible flow within the 
Aerothermopress or . 

:for Professor A # H .Shapiro, Department of Mechanical Engineering, 
DIC 5-6985, GNR I5ori-07&78 

:by B.D.Gavril(DIC), 10 hours 

sDCL : WWI, 111 minutes 

Computations relative to the thermodynamic performance character- 
istics of the aerothermopressor were resumed during the current biweekly 
period* These computations were concerned , with improving performance by 
suitable cross-sectional area, variation and were carried out in a routine 
manner with the existing aerothermopressor WWI program. 

An addition to the program was made which enabled the entropy of 
the flow to be computed and printed with the other properties of the flow. 
Due to the relative length of this subroutine, it was necessary to use the 
.-auxiliary magnetic drum. 

As time permits, it is planned to revise the aerothermopressor 
program while at the same time taking advantage of CS II . ! The current 
program involves the one bank PA and was written in April, 1953* 

.122 B. Coulomb Wave Functions . Regular and irregular solutions of the 

radial Schrodinger equation with, a coulomb potential are being 
sought. These solutions iraist have the proper asymptotic form 
to correspond to scattered waves from such a potential. The 
solutions are written as, a power series in the radial variable 
and the coefficients are determined by recursion formulae* 

:for Professor H.Feshbach, Professor P.M.Morse, Physics Department 

;by A. Terakin( Res .Assist. CMMG), 25 hours 

sDCL: Staff j .'.§iegeli"3 hours; WWI, 19 minutes 

A test program has been written and successfully tested to solve 
the recursion relation: 
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(n-Ul)(#L) a L B « 2 n 4 - a^ 

1 
i 

where L and rj are parameters and n( which starts at -L) is to run until the 
coefficients ar get less than a preas signed amount. These coefficients will 
be eventually needed in the power series expansion of the Coulomb wave 
functions* 

123 C. Earth Resistivity measurements are used to calculate the Slichter 
kernel function which, in special .1 cases, can be analyzed to 
give the actual distribution of resistivity. The method in- 
volves least-square fitting a set of polynomials to the measured 
surface-potential function and integrating the product of this 
set and the zero-order Bessel functipn. 

:for P. M .Hurley, Department of Geology and Geophysics, D1C 5-6915 

:by K, Vozoff(Res.Assist.CMMC), 50 hours 

:DCL Staff: Demur jian, 1 hour; WWI, 1 21 minutes 

The program for calculating the Slichter kernel integral by iter- 
ation was completed, but the integration took far too much time. A shorter 
method of doing it was suggested by T.R.Madden, and this was analyzed and 
programmed. The new method is similar to the old one in that it first fits 
a set of Legendre polynomials to the original potential data. The integral 
is then of the form 

k(M «» A°r v P(r) J (A r) dr +h[rf%fj Q (hr)dr 
«7Q 5 (p) J ( >> r)dr +* [r^ /j ( * r)dr 



where x s= Ar 
a = h r c 



c 

X 

c 



k( h ) =y C C o +C l x 4C 2 x2+ ^^ 4C 3* 5 J J (x)dx + K /j o (x)dx 

The integrals 
ft, 

x <T (x)dx 

can be evaluated analytically. For n odd, the integral appears in closed 
form. For n even, the integral contains the sum of an infinite number of 
terms 

W a) + "W a) + --— + Vi+2i (a) + -] 
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In the range of values of a considered, this sum converges so that terms 
beyond J or/ 3 ) can ^ e neglected. 

The remainder term, 

K J J (x) dx 
o 

is evaluated in the same manner(n - 0) . 

The resulting expression is a sum of the form 
30 
kU) = 2 d.J.(a) 
i=0 x x 

to be evaluated for each value of fN and the corresponding value of a. The 
values J (a) and J,(a) are calculated using the first 12 terms in the series 
expansion 1 , which gives an accuracy of 10~-5 in this range of a. The following 
terms J 2 , J~, etc., are calculated using the recursion expression for theiru 

One error in the program has been found and corrected. Further 
tests have been delayed by conversion difficulties. 

The first half of the program for the analysis of the kernel function 
has been tested. An error has been found at the end of this program and has 
been corrected. 

126 C. A Data Reduction Program for use in the Servo-mechanisms Laboratory 
is being developed in separate stages to be combined at a later 
date. "he first stage is concerned with devising a program to 
fit polynomials to arbitrary empirical functions using a least 
squared error criterion. The procedure makes use of Legendre 
polynomials and matrix multiplication. 
: for J.E.Ward, Servo Laboratory, DIG No. 7138, AF33( 6l6) 2038 
:by P. T. Ross (DIC), 5 hours; Turyn, 40 hours; Hamilton, 60 hours 
:DCL: "WWI, 66 minutes 

Tests run on the Polynomial Fit program seem to indicate that 
smoothness of differences of data does not affect the degree of the opti- 
mum approximating polynomial. For example, the optijjtum polynomial for a 
section of the sine function was of degree 11; it was not changed very 
much when the sine put in was accurate to only 5 digits instead of 8. 

It was discovered that poor runs of part 1 of the Data Reduction 
problem were due to a numerical bias in the data which had not been removed. 
The program has been changed to correct this# 

The scales on some of the scope plots in the Polynomial Fit program 
are being changed* 

132 C. Subroutines for the Numerically Controlled Milling Machine are being 
revised and tested. The set of subroutines facilitates program- 
ming of the computations involved in the preparation of numerical 
data used to control the milling machine. The subroutines 
involve routine numerical and logical operations. 

sfor J.O.McDonough.Servo Laboratory, DIC No. 6873 

:by J. H,Runyon( Res. Assist. E.E.), 20 hours 

:DCL: WI, 57 minutes 
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A subroutine for determining cut spacing across a particular cross- 
section for constant scallop height was successfully tested. 

A series of milling machine tapes for a conical surface is being 
prepared. These illustrate two different methods of cut spacing and two 
degrees of approximation to the surface. If the results of this invest- 
igation are satisfactory, tapes will be prepared for a conic wing section, 

1^9 C ♦ Digital Methods of Detecting Signal from Noise are being investi- 
gated, A sequence of binary numbers will simulate the message 
wherein regions of high density of ones are signal regions and 
those with low density of ones are noise regions. Various 
methods of detecting the change from one region to another, as 
well as the length and midpoint of the signal regions are being 
studied, 

s for J, V. Harrington, Lincoln Laboratory 

:by G.P.Dinneen, Lincoln Laboratory, 20 hours; Reed, 2 hours 

:DCLs WWI, 106 minutes 

The second phase of this study has been started. A successful run 
was made to determine the distribution of ones in the generated intervals 
corresponding to noise and corresponding to signal. Emphasis during this 
phase will be on the modified success run detector. 

Some difficulty, due to a programmer's error, has been encountered 
in the calculation of the mean and variance of the locations of the center 
of the signal regions. 

152 D, Diffusion in an Oxide-Coated Cathode is a program to calculate the 
effects of combined thermal and electrolytic diffusion that 
occur in an oxide-coated cathode when current is caused to 
flow through the cathode, 

:for W.B.Nottingham, Physics Department, DIC No, 63^5 

sty H.B.Frost(Res, Assist, E.E,), , 5 hour 

:DCL: WWI, 66 minutes 

During the past period two runs were completed successfully. One 
of these runs was the completion of an earlier run, while the other run 
was the calculation of a complete decay and recovery characteristic, with 
variable current, constant field conditions. These runs complete this 
problem as it is now envisioned, although it is possible that one or more 
additional runs may be required if the calculated data prove insufficient 
in range for the necessary experimental comparisons, 

155 B, Synoptic Climatology , A multiple regression formula is used to 

predict temperatures 'from pressure distributions described by 
Tschebycheff polynomials. The matrix of scalar products 
which is used in the calculation of the coefficients of the 
multiple-regression system is being calculated on WWT, 
1 for Professor T,,F,Malone, Meteorology Department 
:by R. Miller(DIC), 40 hours; Friedman, *40 hours 
:DCL Staffs Demurjian, 2.75 hours; Denman, 3 hours; Porter, 6.5 
hours; WWI, 419 minutes 

A cross-product program was successfully tested. This program has 
been used to obtain the cross-products of the orthogonal polynomial coeffic- 
ients with the temperature at five localities and the. precipitation ever 
three river drainage (continued on page 9) 
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basins* .Four sets of orthogonal polynomial coefficients were used denoting, 
respectively, the sea-level pressure surface, and the upper-air pressure 
surface using height, probability of height, and standardized height units, 

Friedman 

Specification of the surface pressure pattern for ha If -a -hemisphere 
has been completed for the months of January 19*^8-1952. ^he sums and cross- 
products necessary in a 52 x 52 matrix used to predict surface pressure 
values for a major portion of the North American Continent twenty-four hours 
in advance have been completed. Solution of the matrix is now being performed 
by desk calculator using the Crout method. 

Miller 

161 C. Response of Mass-Plasti c Spring System to Transient Loading ; 2nd 
order non-linear difference equation representing the response 
of building foundations to transient shock loads is being 
studied in order to develop criteria for the design of blast 
resistant foundations. The footing is represented by a concen- 
trated mass, and the soil by a variable mass and an elastic- 
plastic spring, A Runge-Kutta fourth order integration pro- 
cedure will be used, 

sfor Professor R, Whitman, Asst. Prof, of Soil Mechanics 

s by S. Sydney(Res.Asst,CMMC)y 25 hours 

sDCL Staff; Kopley, 1 hour; WC, J6 minutes 

This program, prepared for an analysis of cohesionless soils, such 
as sand and gravel, is operating satisfactorily. Some results have been 
obtained and have been checked by hand calculations. Production runs will 
be made and analyzed for this case before the work on cohesive soils such 
as clays is started, 

163 C, Ferrite Phase Shifters in Rectangular Wave Guide; transcendental 
equation. The electromagnetic boundary value problem dealing 
with the non-reciprocal ferrite phase shifter in rectangular 
wave guide has been solved. Special cases (assuming negligible 
magnetic losses) of the resulting complicated transcendental 
equation have been computed by hand. Additional computations 
by machine are required to investigate other ferrite materials 
and to establish a frequency dependence. Since magnetic loss 
is a figure of merit for the system described it is essential 
that some investigation be made. For cases in which the losses 
are significant the system will be a non-reciprocal ferrite 
attenuator. The numerical solution will be obtained by oper- 
ating on two simultaneous transcendental equations, 

sfor Dr* Benjamin Lax (DDL) 

sby K.J.Button(DDL), 60 hours 

sDCL Staffs Demur jian, 1 hour; WWI, 28 minutes 

The fourth routine for the first part of the problem was run using 
parameters corresponding to a new ferrite and a different frequency that 
is now being used by the experimenters, A plot of the results revealed a 
small discontinuity of the first kind (finite jump) after about 2/3 of the 
program had been performed. When a detailed analysis failed to detect any 
error in the program, a fifth run was attempted. This fifth performance 
ran correctly and did not reproduce the discontinuity. A tape was prepared 
which doubled the thickness of the -ferrite and a sixth performance was requested, 
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A second test performance .was submitted for the second part of the 
problem, ^he program parameters printed out for purposes of analysis were 
hand calculated and found to be correct. The end results were of the proper 
order of magnitude. Since the test program was judged to be correct, addit- 
ional routines were added to increase the accuracy and quantity of the results 
The third performance printed results that appear to be of the correct order 
of magnitude. This performance utilized parameters corresponding to the 
new ferrite and frequency described in Part I. 

Future Plans: Since the sixth performance of Part I (just received) 
has been found to be correct, additional tapes will be made to change para- 
meters for the investigation of resonance conditions, frequency band widths, 
new ferrites and new frequency bands. The same conditions will be imposed 
upon Part II after a rigorous verification of the accuracy of the results 
has been calculated, 

166 C. Construction and Testing of a Delta Wing Flutter Model is being 
effected by replacing the actual wing by a structurally 
equivalent lattice network. An iterative procedure involving 
the evaluation of a matrix equation has been evolved for deter- 
mining the bending and torsional stiffnesses of the component 
members of the network. 

jfor M.M.Chen(DIC) 

:by S.Gravitz( Res .Assist. Aero.Eng»), 40 hours 

sDCL Staff: Porter, 3.5 hours; WWI, 18 minutes 

The routine for formulating a matrix with each element a sum of 
products has been completed and successfully tested. 

Routines for adding matrices with each element a sum of coefficients 
of unknown terms and for equating two matrices term-by-term to form a set 
of simultaneous algebraic equations were observed to entail minor modific- 
ations of subroutines previously developed for this problem. 

A routine for matrix multiplication with one matrix having elements 
composed of coefficients of unknown quantities has been completed and is 
in the process of being tested. An extension of this routine to matrices 
too large to fit entirely into high speed storage is being evolved. 

The program for the future will consist of coordinating the sub- 
routines previously evolved into an integrated. program for the ultimate 
purpose of solving the wing simulation problem. 

168 D. Indicial Downwash behind a Two-Dimens ional Wing . In the analysis 
of the response of an airplane to a sharp-edge gust, and par- 
ticularly in the calculation of stresses in the horizontal 
tail, it is important to know the downwash at the tail caused 
by the lift response of the wing to the gust. In the present 
solution for the downwash behind a two-dimensional wing, some- 
thing more than the "indicial downwash" is sought. This is 
, effected by allowing the gust front itself to have an arbitrary 
horizontal velocity which, combined with the velocity of the 
airplane, U, results in the wing's penetrating the gust front 
at a velocity V. 

:for Professor, H. Ashley, Aero. Eng. Dept. , DIC 6727 

:by N.P.Hobbs, 60 hours 

:DCL Staff: Porter, .5 hour; WWI, 18 minutes 



M-27^2 page 12 

The function Cj(A) (see Biweekly Report of March 7. 195*0 has been 
computed on Whirlwind and found to be unity"* The majority of the past two 
weeks has accordingly been spent in proving: rigorously that Gj(k) fc 1« 
Having proved this to be true, the outer integral involving Cj(A) was eval- 
uated analytically. 

The remaining problem is to evaluate numerically Cjj(k t X*) and the 
outer integral involving Cjj(A,X*)» 

171 G. Improved Power Spectrum Estimates are to be obtained by investigating 
modifications to existing techniques • for the numerical calculation 
of Fourier cosine transforms to minimize the effects of truncation 
and to supply a confidence curve, based upon the internal varia- 
tions of the calculations, which will help to evaluate the signif- 
icance of the resulting spectrum estimate. The object is to 
obtain a method which will give significantly better power spectrum 
estimates than can be obtained using existing techniques. 

s for J. E.Ward, Project Engineer, Servome onanisms Laboratory 

:by D,T. Ross (DIC), 100 hours; Hamilton, 20 hours 

sDCLs WWI, 126 minutes 

In the proposed procedure, a series of spectrum estimates, called 
the First Spectrum, Second Spectrum, etc., will be calculated. Each spectrum 
is closer to the true spectrum than the previous one by a factor which is 
proportional to the truncation length, T. The unnormalized first spectrum 
S-. (w,T) is the standard estimate without any modification 

S. (w,T) s 2 J f(t)cos wtdt, 

° th 

where f(t) is the function to be transformed. The unnormalized n Spectrum 

is recursively defined in terms of the (n-l)££ Spectrum by 

S n *(w,T) 'J tS , n * 1 (w,t)dt. 
The normalized spectra are then defined by 

S n (w,T) =^f^-l S n *(w,T). 

Thus it is theoretically possible, for a given truncation length T, to obtain 
a spectrum estimate for the given section of the function (from to T), -sihich 
is accurate to any specified tolerance, merely by taking n large enough. 

The measure of indecision which will be tested as part of this study 
will be the amount of oscillation in the process over and above the oscillation 
which is normal for the function cos\)t. Possible forms for this measure 
of indecision which are easy to compute at the same time as the spectra are 
computed are 

,T 



^ / |cos 2 vt - f(t) cost>t| dt 



and 
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- / (cos 2 yt-f(t) cos ut) dt 

where, again, f(t) is the function being transformed* Both forms will be 
tested* 

Every method for obtaining spectra by computational means can give only 
estimates of the true spectrum. Hence, such a process may be considered 
to have a certain filter characteristic , which we call the M spectral 
window" of the process* The spectral window of any such process when the 
frequency V . is being investigated may be taken as 

T 

u (w) = 2- \S^ ) cos kV cos kw 
o 

where k is the index of the digital process (corresponding to the time var- 
iable t) and hjj(v ) includes all coefficients and special operators of the 
process which are used when the frequency V is the scanning frequency. Thus 
u^ (w) measures not only the convolving function (e.g.* (sin x)/x) for 
ordinary truncation, but also the effects of the actual numerical processes 
used. Note that it is obtained by inserting the function cos \)t into the 
process and hence is easily obtained. All stages of this investigation 
will be tested by finding u v (w) but final evaluations will be based upon 
the performance with complicated sums of cosines, which approximate auto- 
correlation functions actually found from nature. 

For testing purposes, a program to generate functions of the form 

YL 

2- A(w.)cos w.t 

for an arbitrary number of frequencies, w^, has been written. The present 
version allows the use of up to fifty frequencies with arbitrary amplitudes. 
It will thus be possible to approximate the types of autocorrelation 
functions actually found from natural processes by functions with known 
spectra . 

The improved Fourier Transform program using weighted means appears to 
be operative. Several minor logical and technical errors were found and 
corrected, ^he second (equivalent) program which uses pre-multiplicatiOn 
by suitable weighting functions must still be tested. The program to generate 
sums of cosines for use as test functions has been tested and the latest 
version should be correct. Kow several test runs will be made to determine 
the characteristics of the programs* the method and to establish sound 
procedures for general application. 

172 Bk Overlap Integrals of Molecular and Crystal Physics . Two-center 
overlap integrals are to be evaluated between various Slater 
atomic orbitals, which are of the form: (power of r) x 
(exp -ar) x (spherical harmonic) . By use of prolate spher- 
oidal coordinates, formulas in terms of simple functions can 
be derived for these integrals but are unsatisfactory because 
they are of formidable complexity and have false singularities* 
The integrals are evaluated by recasting the expressions in 
terms of spherical Bessel functions of imaginary argument, 
which can be generated by a high-speed computer. 
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:for Professor J ,C. Slater, Physics Department 
:by F.J.CorbatoXRes.Asst.CMMC), 50 hours 
:DCL; WWI, 90 minutes 

In the Tight-Binding approximation applied to a solid, two-center 
integrals are the most important quantities to be evaluated. Heretofore, 
these evaluations have been a major limitation of the method. Although 
in principle, closed analytic forms have been developed for these integrals, 
(C.C.J.Roothaan, Jour, of Chem. Phys., 12., 1445,(1951)) in practice the 
formulas are not only formidable but are very poor for computational pur- 
poses because of subtractions of large numbers and the resultant low 
accuracy* These integrals have been the object of much analytical work in 
the past, and if capable of being evaluated in a simple fashion, would be 
of great interest to physicists especially in the fields of S olid-State and 
Molecular Structure. 

The basic subroutines of the problem have been coded and tested 
using CS II* The final program for the overlap integral subroutine has 
been written but must still be typed and given a final check. The more 
educational programming mistakes were: 

1) Any cycle instruction with a buffer address will not operate 
correctly. 

2) The subroutine LSR FU1 for e*, calculates e^ x * regardless of the 
sign of x, and then divides the result into 1 if x is negative* Consequently, 
for | x \ ^ about +50 • , even with x negative , the routine fails on an 
interpretive overflow during an its2t order. One solution of this difficulty 
is to modify the subroutine by replacing the 2t addresses by a buffer address* 
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1.3 Operating Statistics 
1«31 Computer Time 

The following indicates the distribution of WI time allocated to the 
S&EC Group. 



Programs 
Conversion 
Magnetic Drum Test 
Magnetic Tape Test 
Scope Calibration 
Demonstrations(#13l) 



70 hours, 23 minutes 
2k hours, 06 minutes 
25 minutes 
kl minutes 
Itfy minutes 
Ik minutes 



Total Time Used 
Total Time Assigned 
Usable Time, Percentage 
Number of Programs 

1.32 Program Time Distribution 



96 hours, 29 minutes 
216. hours, JO minutes 

3Jk 



The following table attempts to show how the ffiH. time expended on S&EC 
programs was distributed with respect to machine runs that gave meaningful 
results (productive computer time) and runs that gave unsatisfactory results 
(lost computer time). Productive computer time is subdivided to indicate 
the time involved in actual computations as contrasted with the time expended 
getting information out of WWI. Computer time lost is subdivided to show 
the portion of time lost, due to errors in the programmers formulation of 
his problem (logical errors); due to errors in the programmer's use of the 
WE code, . CS. Conventions, etc. (technical errors) ; due to tape preparation 
errors?; due to errors by the S&EC computer operators in running the program; 
due to malfunctioning of terminal equipment; and finally, due to miscellan- 
eous causes. 

These times are determined as percentages of the time listed above in 
section 1.31 for programs. The times used in computing these figures are 
extracted from the biweekly report forms submitted by the various program- 
mers who have used S&EC allocated WWI time. 

1. Productive Computer Time 
Computation 48.1$ 
Output 12.7$ 

2. Computer Time Lost Due to Programmers Errors 
Technical 17*3$ 
Logical 12. 0$ 

3. Computer Time Lost Due to Other Difficulties 
Tape Preparation 2.0$ 
Operators' Errors 1.*$. 
Terminal Equipment Malfunction 4.0$ 
Miscellaneous 2.5$ 

1*33 Tape Preparation 

An aVtempt is being made to obtain some idea of. the time expended in 
the preparation of tapes. During the past biweekly period a check was made 
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on the tapes processed. 

Due to the variations in procedures involved we have distinguished 
among original complete tapes and the following three types: typed mod- 
ifications - changes of 11 or more registers which must be ,typed, conver- 
ted and then attached to the main program or changes which must be made 
in the body of a Flexowriter tape; manual modifications - changes punched 
directly in 556 form and attached to a converted tape; combined tapes - 
which require duplication of two or more complete tapes* 

The following information was compiled: 

Complete Typed Manual Combined 
Tapes Mods Mods Tapes 



No. of Tapes 103 67 46 14 

No. of Registers 21606 6048 340 

Time Consumed 58hrs.l3min. 51hrs.59min. 7hrs. 45min. 6hrs. 3min. 

Thus, it may be seen that the average length of an original complete 
tape is 167.7 registers requiring 33.9 minutes to prepare. A typed 
modification averages 13.2 registers in length and requires 25.9 minutes 
to prepare which Manual Modifications average 7,4 registers and require 
10.1 minutes for preparation. 

In addition, to the time accounted for above, a total of 12 hours and 45 , 
minutes was spent printing computer' results for photographic reproduction. 

2. ACADEMIC PROGRAM 

CS Programming Course 

The CS introductory two-week programming course began on March 15 
with an enrollment of 18 students. These students represent the following 
groups: Aeroelastic and Structures Research Laboratory, Instrumentation 
Laboratory, MIT Electrical Engineering Department, MIT Chemical Engineering 
Department, Dynamic Analysis and Control Laboratory, MIT Mathematics 
Department, Lincoln Laboratory, Geophysical Analysis Group, and the Lab- 
oratory for Nuclear Science and Engineering. 

Several chapters of the Manual being used for the first week's 
work are being rewritten to conform with CS II conventions. It is expect- 
ed that this manual will have been revised before the course is given again. 

Seminar on Advanced Programming Techniques of the Digital Computer 
Laboratory Staff . 

On Friday, March 12, Dr. H. H. Denman discussed the uses of magnetic 
tape memory and its associated delayed printer equipment. 

Mr. P.R.Bagley described the Memory Test Computer on March 19. 

On Friday, April 9# Professor C.¥.Adams will talk on "Recent Devel- 
opments in the Digital Computer Field Including Current Trends and New 
Computers" « 
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S&EC Movie 

The script for the movie "Making Electrons Count" has been rewritten 
and "shooting" will begin on the revised version of the film within a few 
days • 

Industrial Liaison Conference 

A conference on "Control Applications in Business and Industrial 
Systems" will be held April 5-6" at MIT under the sponsorship of MIT's 
Industrial Liaison Office. Those who plant to attend should contact 
Miss Wellington at extension 2693 immediately. 

Seminar on Computing Machine Methods 

On Tuesday, March 9» E.A.Robinson, director of the Geophysical 
Analysis Group at MIT, spoke on "Petroleum Exploration and Digital 
Computers. 

On Tuesday, April 6, Professor James B. Reswick of the MIT Depart- 
ment of Mechanical Engineering will speak on "An Electronic Analog Device 
For Delay Line Synthesis" • The seminar is held at k FM in Room 12-182. 
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3» COMPUTER ENGINEERING 

3.1 WWI System Operation (L.L .Holmes, A. J .Roberts) 

A prototype of the new core-memory sense amplifier is nearly complete. 
It is planned to install the unit in WWI for tests on Tuesday, 16 March. 
If the results of the tests are satisfactory, the unit will be left in ser- 
vice for a reliability check* 

All of the fixed- voltage switching panels in Rack PO have been replaced 
by a new model. This is the initial step in our plan to replace all of the 
fixed-r voltage switching panels and voltage- variation switching panels in 
P Row. The models contain more reliable relays and should therefore ease 
our relay maintenance work load. 

During the past biweekly period, the majority of the systems testing 
time was used to investigate the main control section of the computer in 
search of the origin of several transient transfer-check alarms that have 
occurred. Several faulty tubes and crystals were located, but the trouble 
is believed to be still present. More testing will be done during the 
systems-testing period. 

3.11 Typewriters and Paper Tape (L.H.Norcott) 

Commercial Controls Corp. reports that our three new Flexcwriters are 
now going through their final assembly line and will be shipped to us short- 
ly. As soon as they are received, the FXexo shop will modify them for use 
with Whirlwind. 

3.2 Terminal Equipment (R.H.Gould) 

One K1084P?M 16-inch cathode-ray tube that was removed from service 
because of spurious eB&ssion from the control grid has been successfully 
rehabilitated. R-*F bombing of the control grid stopped the grid emission 
without apparent damage to the cathode. 

Logical difficulties with block control have recently come to light, 
and a rather basic redesign of block control may be desirable. Other 
changes to in-out control to simplify magnetic-tape operation and to simplify 
in-out control generally are being planned. 

Work continues on bringing in-out drawings up to date. 

3.21 Buffer Drum System (K.E.Mc Vicar) 

The buffer drum is now connected to the computer at all times and is 
available to programmers without special arrangement. 

The auxiliary-storage section of the buffer drum does not yet have a 
parity digit, but it is. checked daily and can be used with reasonable 
assurance of reliable operation. The interlace on the auxiliary-storage 
section of the buffer drum is eight instead of four as on the auxiliary 
drum. It takes about twice as long to make a complete block transfer of 
one group from the buffer as it does from the auxiliary. 
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3»22 Magnetic Drums (H.L.Ziegler) 

Installation of electronic-write switching of heads in the auxiliary 
drum continues to progress satisfactorily. All hardware is on hand, and 
most of it is installed. Filament and d-c power wiring are nearly complete, 
and the bulk of the signal leads have been made up into cable form. These 
leads are to be placed in the wiring ducts during installation time Monday, 
15 March, after which final wiring of the chassis plugs can begin almost 
immediately. 

The magnetic- drum-PETR monitor system is also nearing completion. As 
presently planned, monitoring will be available for the entire auxiliary 
drum, groups 4-7 of the buffer drum, and the PETR output channels. Included 
in the new system is an entirely new switching panel designed to consolidate 
the present widely scattered controls and also to provide some desirable 
interlocking of these controls. 

3 «23 Magnetic Tape (E. P. Farnsworth) 

The unit 2 delayed print-out equipment and test programs have been check- 
ed out and can be plugged into the second Flexowr iter when the enclosure 
is completed. Switches to interchange unit 2 and unit 3 IOC, unit 2 and 
unit 3 print-outs, and to convert the unit 2 Flexowriter to standard FL 
or delayed print-out are under construction in the shop. The rewind indi- 
cator gong is installed and operating, and "ON LIMED" lights are being 
mounted on each tape-unit control panel. Modifications are being made to 
improve the two-tone chime effect to prevent confusion with the computer 
alarm, and to change the "Rewind Unit M switch near PE3K from a toggle to 
a push-button switch. 

The delayed print-out functional schematic drawing was completed, and 
the block schematic for the entire magnetic- tape system has been brought 
up to date. Both drawings are in the WWI service file. 

3«24 Flexowriter Equipment (M. Demur jian) 



Discussions were held with Norcott and Garrill on the grounding of 
Flexowriter equipment,. It was decided that all Flexowriters would have 
special Hubbell plugs attached and the wall receptacles would be replaced 
by Hubbell type wherever Flexowriters are used. For the mobile uses 
adapters would be available. - 

3 • 25 Magnetic Tape Equipment (M. Demur jian) 

The new Delayed Output box is in Test Control. As soon as the control 
boxes and cabies r are completed, Al Perry will connect it to the Magnetic 
Tape Units. The expected date of delivery for these is April 2, 1954. 

Perry located and remedied the cause of spurious digits and characters 
which were appearing in the delayed output. ' 
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h * ADMINISTRATION AND PERSONNEL 

New Non-Staff (R. A. Osborne) 

Renee Feinstein is a new secretary in the Publications Office. 
Lester Gediman is a new member of the Drafting Department . 

Terminated Non-Staff 

Jean Pfaff 
Mildred Stickney 
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